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We present the elements of QED®QCD exponentiation and its interplay 
with shower/ME matching in precision LHC physics scenarios. Applications 
to single heavy gauge boson production at hadron colliders are illustrated. 
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In the LHC environment, precision predictions for the effects of multiple 
gluon and multiple photon radiative processes will be needed to realize the 
true potential of the attendant physics program. For example, while the cur- 
rent precision tag for the luminosity at FNAL is at the ~ 7% level [1], the 
high precision requirements for the LHC dictate an experimental precision 
tag for the luminosity at the 2% level [2]. This means that the theoretical 
precision tag requirement for the corresponding luminosity processes, such 
as single W,Z production with the subsequent decay into light lepton pairs, 
must be at the 1% level in order not to spoil the over-all precision of the 
respective luminosity determinations at the LHC. This theoretical precision 
tag means that multiple gluon and multiple photon radiative effects in the 
latter processes must be controlled to the stated precision. With this objec- 
tive in mind, we have developed the theory of QED ® QCD exponentiation 
to allow the simultaneous resummation of the multiple gluon and multiple 
photon radiative effects in LHC physics processes, to be realized ultimately 
by MC methods on an event-by-event basis in the presence of parton showers 
in a framework which allows us to systematically improve the accuracy of the 
calculations without double-counting of effects in principle to all orders in 
both a s and a. 

Specifically, the new QED ®QCD exponentiation theory is an extension 
of the QCD exponentiation theory presented in Refs. [3] 1 . We recall that in 
the latter references it has been established that the following result holds 
for a process such as q + q' — > V + n(G) + X — > ££' + n(g) + X: 



where gluon residuals j3 n (ki, . . . , k n ) , defined by Ref. [3], are free of all in- 
frared divergences to all orders in a s (Q). The functions SU Mm(QCD), Dqcb, 
together with the basic infrared functions Bq^ d , Bq^ d , Sq 1 q D are specified 
in Ref. [3]. Here V = W ± ,Z, and £ = e,fi, £' — u e , z^(e, /x) respectively for 

1 We stress that the formal proof of exponentiation in non-Abclian gauge theories in the 
eikonal approximation is given in Ref. [4] . The results in Ref. [3] are in contrast exact but 
have an exponent that only contains the leading contribution of the exponent in Ref. [4] . 
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V = W + (Z), and £ = v e ,v u , £' = e,/i respectively for V = W~. We call 
attention to the essential compensation between the left over genuine non- 
Abelian IR virtual and real singularities between J dPh(3 n and J dPh(3 n+ \ 

respectively that really allows us to isolate (3j and distinguishes QCD from 
QED, where no such compensation occurs. The result in (P) has been real- 
ized by Monte Carlo methods [3]. See also Refs. [5-7] for exact 0(a 2 s ) and 
Refs. [8-10] for exact 0(a) results on the W,Z production processes which 
we discuss here. 

The new QED ® QCD theory is obtained by simultaneously exponenti- 
ating the large IR terms in QCD and the exact IR divergent terms in QED, 
so that we arrive at the new result 

d& m = e SUM IR (QCED) 
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where the new YFS [11,12] residuals, defined in Ref. [13], (3 n>m (ki, . . . , k n ; k[, . . . , k' r 
with n hard gluons and m hard photons, represent the successive appli- 
cation of the YFS expansion first for QCD and subsequently for QED. 
The functions SUMir(QCED), Dqced are determined from their analogs 
SUMir(QCD), Dqcb via the substitutions 

Tjnls Tjnls I Tjnls jjnls 

n QCD ~ y D QCD "+~ n QED = n QCED> 

fynls finis , finis finis 

n QCD ~^ n QCD ' n QED = D QCEDi 

nnls nnls , nnls nnls {q\ 

S QCD ~^ S QCD ' a QED = ^QCED \ 6 ) 

everywhere in expressions for the latter functions given in Refs. [3]. The 
residuals /? n , m (/ci, . . . ,k n ;k[, . . . , k' m ) are free of all infrared singularities and 
the result in (J2J) is a representation that is exact and that can therefore be 
used to make contact with parton shower MC's without double counting 
or the unnecessary averaging of effects such as the gluon azimuthal angular 
distribution relative to its parent's momentum direction. 

In the respective infrared algebra (QCED) in (J2J), the average Bjorken x 
values 

x avg (QED) = j(QED)/ (1 + >y(QED)) 
x av9 (QCD) = j(QCD)/(1 + -/(QCD)) 



mil 



3 



where 7(A) = ^f±(L s - 1), A = QED, QCD, with C A = Q 2 f ,C F , respec- 
tively, for A = QED, QCD and the big log L s , imply that QCD dominant 
corrections happen an order of magnitude earlier than those for QED. This 
means that the leading /?o,o-level gives already a good estimate of the size of 
the interplay between the higher order QED and QCD effects which we will 
use to illustrate (0) here. 

More precisely, for the processes pp —> V+n(j)+m(g)+X — > ££'+n'('y) + 
m(g) + X, where V = W ± , Z,and £ = e,/i, I' — v e , i/ M (e,//) respectively for 
V = W + (Z), and £ = v e , v^, £' = e, \i respectively for V = W~ , we have the 
usual formula (we use the standard notation here [13]) 



and we use the result in (J2J here with semi-analytical methods and structure 
functions from Ref. [14]. A Monte Carlo realization will appear elsewhere [15]. 

We do not attempt in the present discussion to replace HERWIG [16] 
and/or PYTHIA [17] - we intend here to combine our exact YFS calculus 
with HERWIG and/or PYTHIA by using the latter to generate a parton 
shower starting from the initial (xi,X2) point at factorization scale fi after 
this point is provided by the {Fi}. This combination of theoretical constructs 
can be systematically improved with exact results order- by-order in a s , where 
currently the state of the art in such a calculation is the work in Refs. [18] 
which accomplishes the combination of an exact 0(a s ) correction with HER- 
WIG. We note that, even in this latter result, the gluon azimuthal angle is 
averaged in the combination. We note that the recent alternative parton 
distribution function evolution MC algorithm in Refs. [19] can also be used 
in our theoretical construction here. Due to its lack of the appropriate color 
coherence [20], we do not consider ISA JET [21] here. 

To illustrate how the combination with Pythia/Herwig can proceed, we 
note that, for example, if we use a quark mass m q as our collinear limit regu- 
lator, DGLAP [22] evolution of the structure functions allows us to factorize 
all the terms that involve powers of the big log L c = ln/z 2 /m 2 — 1 in such 
a way that the evolved structure function contains the effects of summing 
the leading big logs L = In fi 2 / /1q where we have in mind that the evolution 
involves initial data at the scale fiQ. The result is therefore independent of 
m q for m q [ 0. In the context of the DGLAP theory, the factorization scale 
ji represents the largest p±_ of the gluon emission included in the structure 



da exp ( PP ^V + X^U' + X') 
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function. In practice, when we use these structure functions with an exact 
result for the residuals in J2J) , it means that we must in the residuals omit the 
contributions from gluon radiation at scales below fi. This can be shown to 
amount in most cases to replacing L s = Ins/ m 2 q — 1 — > L n \ s = In s/fi 2 but in 
any case it is immediate how to limit the pt in the gluon emission 2 so that 
we do not double count effects. In other words, we apply the standard QCD 
factorization of mass singularities to the cross section in (j2J) in the standard 
way. We may do it with either the mass regulator for the collinear singu- 
larities or with dimensional regularization of such singularities - the final 
result should be independent of this regulator. This would in practice mean 
the following: We first make an event with the formula in which would 
produce an initial beam state at (xi, X2) for the two hard interacting partons 
at the factorization scale \i from the structure functions {Fj} and a corre- 
sponding final state X from the exponentiated cross section in da exp (xiXjs) ; 
the standard Les Houches procedure [23] of showering this event (xi,X2,X) 
would then be used, employing backward evolution of the initial partons. 
If we restrict the pr as we have indicated above, there would be no double 
counting of effects. Let us call this px matching of the shower from the back- 
ward evolution and the matrix elements in the QCED exponentiated cross 
section. 

However, one could ask if it is possible to be more accurate in the use of 
the exact result in (J2J)? Indeed, it is. Just as the residuals j3 n>m (ki, . . . , k n ; k[, . . . , k' m )aie 
computed order by order in perturbation theory from the corresponding ex- 
act perturbative results by expanding the exponents in (J2J) and comparing 
the appropriate corresponding coefficients of the respective powers of a n a™, 
so too can the shower formula which is used to generate the backward evo- 
lution be expanded so that the product of the shower formula's perturba- 
tive expansion, the perturbative expansion of the exponents in (J2J), and the 
perturbative expansions of the residuals can be written as an over-all ex- 
pansion in powers of a n a™ and required to match the respective calculated 
exact result for given order. In this way, new shower subtracted residuals, 

{Pn,m(ki, . . . , k n ; k[, . . . , k' m )}, are calculated that can be used for the entire 
gluon pt phase space with an accuracy of the cross section that should in 
principle be improved compared with the first procedure for shower matching 
presented above. Both approaches are under investigation. 

Returning to the general discussion, we compute, with and without QED, 

2 Hcrc, wc refer to both on-shell and off-shell emitted gluons. 
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the ratio r exp = a exp / 'asorn, where we do not use the narrow resonance ap- 
proximation; for, we wish to set a paradigm for precision heavy vector boson 
studies. The formula which we use for a Born is obtained from that in (j3J) by 
substituting da Born for da cxp therein, where da Born is the respective parton- 
level Born cross section. Specifically, we have from (JTJ) the /^co-level result 



r^max 

a exp (xix 2 s) = / dv"i Q cEDV lQCED ~ l F Y Fs(lQCED)e 8YFS : a B orn(0--v)xix 2 s) 
Jo 

(5) 

where we intend the well-known results for the respective parton-level Born 
cross sections and the value of v max implied by the experimental cuts under 
study. What is new here is the value for the QEDcg>QCD exponent 

Iqced = {2Q}^ + 2CV^} L nls (6) 

where L n [ s = \nxiX2s/ fi 2 when \x is the factorization scale. The functions 
Fyfs{1qced) and S Y fs(1qced) are well-known [12] as well: 



p-lQCEDlE 

Fyfs(iqced) 



r(l + Iqced) ' 
Syfs(1qced) = -.Iqced + (Qf~ + C F — )(2C(2) - -) 
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where ((2) is Riemann's zeta function of argument 2, i.e., 7r 2 /6, and 7^ is 
Euler's constant, i.e., 0.5772... . Using these formulas in (j^EJ) allows us to get 
the results 

(1.1901 , QCED = QCD+QED, LHC 

_ _ J 1.1872 ,QCD, LHC (H] 

rex P ~ S 1.1911 , QCED = QCD+QED, Tevatron W 

k 1.1879 ,QCD, Tevatron. 

We see that QED is at the level of .3% at both LHC and FNAL. This is stable 
under scale variations [13]. We agree with the results in Refs. [5,6,8-10] on 
both of the respective sizes of the QED and QCD effects. The QED effect is 
similar in size to structure function results found in Refs. [24-28], for further 
reference. 

We have shown that YFS theory (EEX and CEEX) extends to non- 
Abelian gauge theory and allows simultaneous exponentiation of QED and 
QCD, QED®QCD exponentiation. For QED®QCD we find that full MC 
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event generator realization is possible in a way that combines our calculus 
with Herwig and Pythia in principle. Semi-analytical results for QED (and 
QCD) threshold effects agree with literature on Z production. As QED is at 
the .3% level, it is needed for 1% LHC theory predictions. We have demon- 
strated a firm basis for the complete 0(a 2 s , aa s , a 2 ) results needed for the 
FNAL/LHC/RHIC/ILC physics and all of the latter are in progress. 
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